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Bound on induced gravitational wave background 
from primordial black holes 
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The today's energy density of the induced (second order) gravitational wave background in the frequency 
region ~ 1CF 3 — 10 3 Hz is constrained using the existing limits on primordial black hole production in the 
early Universe. It is shown, in particular, that at frequencies near ~ 40 Hz (which is the region explored by 
LIGO detector), the value of the induced part of Qgw cannot exceed (1 — 3) x 10~ 7 . The spread of values of 
the bound is caused by the uncertainty in parameters of the gravitational collapse of black holes. 
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Q ! 1. INTRODUCTION 

It is well known now that gravitational waves (GWs) 
'Can be effectively generated by density perturbations 
^^during the radiation dominated era. Tensor and scalar 
O 'perturbations are decoupled at the first order, but it 
£j is not so in higher orders of cosmological perturbation 
^ 'theory. Namely, the primordial density perturbations 
!— 1 and the associated scalar metric perturbations gener- 
ate a cosmological background of GWs at second order 
through a coupling of modes [TJ |2 |3J. In particular, a 
second order contribution to the tensor mode, hj. , de- 
pends quadratically on the first order scalar metric per- 
turbation, i.e., the observed scalar spectrum sources the 
generation of secondary tensor modes. By other words, 
the stochastic spectrum of second order GWs is induced 
by the first order scalar perturbations. Calculations of 
flew at second order and discussions on perspectives of 
measurements of the second order GWs are contained 
in works 01 El El HEl. 

It is natural to conjecture that the detection of GWs 
from primordial density perturbations on small scales 
(not directly probed by observations) could be used to 
constrain overdensities on these scales. However, at the 
present time, gravitational wave background (GWB) is 
not yet detected. So, on the contrary, one can constrain 
GWB using existing limits on amplitudes of primordial 
density perturbations. Such limits are available, in par- 
ticular, from studies of primordial black hole (PBH) pro- 
duction at the beginning of radiation era. 

It is generally known that PBHs form from the den- 
sity perturbations, induced by quantum vacuum fluctu- 
ations during inflationary expansion. The details of the 
PBH formation had been studied in [9l [10] , the astro- 
physical and cosmological constraints on the PBH den- 
sity had been obtained in many subsequent works (see, 
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e.g., the recent review [H]). The detailed constraints on 
amplitudes of density perturbation spectrum had been 
reviewed in recent works [T2l H3] . 

Having the bounds on the primordial density spec- 
trum, it is possible to determine the corresponding 
bounds on amplitudes of produced GWB. 

The plan of the paper is as follows. In Sec. 2 wc 
present the basic relations connecting the primordial 
power spectrum of density perturbations and induced 
GWB. The main result of the paper - the PBH bound 
on the energy density of GWB - is shown in Fig. [1] In 
Sec. 3 the mini-review of theoretical estimates for flow 
is given. Sec. 4 contains the conclusions. 

2. DERIVATION OF THE PBH BOUND 

The energy density of GWs per logarithmic interval 
of k in units of the critical density is given by 
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Ogw(fc, r) — — 



12 \<i(t)H(t) 



V h (k,T). (1) 



The power spectrum Vh is defined by the standard ex- 
pression 

(Mr)Mr)) = l^-S 3 (k + k')V h (k,r), (2) 

where h^fr) is the Fourier component of the tensor met- 
ric perturbation. For scalar-induced GWs, Vh{k,r) is 
obtained from the formula |6J 



V h (k,T)= J dk J ^-p*(|k-k|)7 5 *(fc)J r (fc,fc,//,r), 

-1 

(3) 

where "P*(fc) is the power spectrum of the Bardeen po- 
tential, and J- is the function determined by the transfer 
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functions of the Bardeen potentials and Green functions 
of the evolution equation for the GW amplitude, 

h£ + 2Wi' k + k 2 h k = S(k, t). (4) 

The source term in this equation is given by 

S(k, t) = / d 3 k fc 2 (l - M 2 ) /(k, k, r)* t _ t * £ , (5) 



where / can be explicitly expressed through the Bardeen 
potential's transfer functions [6j. It is seen from Eqs. 
((3j [5]) that both the source term and the power spec- 
trum of induced background of second order GWs de- 
pend quadratically on the first-order scalar perturbation 

It is also seen from the Eq. {3| that for scalar- 
induced GWs, the single mode in scalar spectrum does 
not correspond to the only one mode in Vh- For exam- 
ple, for the J-function-like spectrum Vn(k) ~ 8(k — ko), 
the GW spectrum turns out to be continuous and 
stretching from to 2kg [5]. However, the order of mag- 
nitude of wave numbers of induced GWs is typically the 
same as of scalar perturbations [SJ, so, the GWs typical 
wave number that will be generated from perturbations 
entering horizon at its mass scale Mh can be estimated 

by 



k — k e 



( M h 



-1/2 



9* 

,9*eq , 

2 n 10 23 (M h [g\y 1/2 Mpc" 



-1/12 



(6) 



where in the last equality we have adopted that g 9eq « 3, 
g* » 100, M eq = 1.3 x 10 49 g • (O m /i 2 )- 2 « 8 x 10 50 g. 
The connection between / and k for GW is 



f=7T- = 1-54 x 10~ 15 ( — — r 
2tt VMpc" 1 



Hz. 



(7) 



If PBHs form from a scalar spectrum of perturba- 
tions at a horizon mass scale Mh, the typical PBH mass 
will be of order of Mh (see, e.g., [H]), so ((BJ relates the 
typical PBH mass with the characteristic k of second 
order GWs produced. 

We assume in this work, that the power spectrum 
Vfiik) can have a peak at some wave number ko. Such 
peaks can arise in several types of inflationary mod- 
els (see, e.g., [14] and references therein) and lead to 
PBH production, so that the peak's parameters (such 
as height and width) can be constrained from PBH 
nonobservation. It is convenient to use some kind of 
parametrization to model the realistic peaked power 
spectrum of finite width, e.g., 



hVn{k) = B + (\gVl-B)exv 



Qgk/ko) 



where B —8.6, V\ and E characterize the height and 
width of the peak, fco is the position of its maximum. 
Parameters of such a distribution have been constrained 
previously [12| from non-observation of PBHs and prod- 
ucts of their Hawking evaporation. 

The frequency region ~ 10 3 — 10 -3 Hz which we 
explore in this work corresponds to horizon and PBH 
mass region of ~ 10 11 — 10 23 g. For a rather wide peak 
(E > 1), all constraints on V\ for this region are in 
the range 0.01 — 0.04, depending on k and the param- 
eter S c , which is a density contrast threshold for PBH 
production (see [T3] for details). 

Particularly, it turns out that for mass range ~ 
10 13 — 10 17 g, the best constraints are obtained from 
the requirement that the diffuse gamma ray flux pro- 
duced by PBHs does not exceed the observed extra- 
galactic one. The condition that the fraction of the en- 
ergy density of the universe contained in PBHs, SIpbhj 
does not exceed the one for non-baryonic dark mat- 
ter, r2 n bm ~ 0-2, gives the most stringent limit for 
Mbh ^ 10 17 g (this constraint exists for black holes 
with initial mass Mbh > M* ~ 5 x 10 14 g, i.e., ones 
that did not evaporate up to the present time). For 
Mbh > 10 13 g we use constraints on obtained in 
[12] using the above two conditions. 

In the region Mbh ~ 10 11 — 10 13 g several types of 
constraints coexist, including ones considering PBH in- 
fluence on the photodissociation of deuterium [T5J [16] , 
CMB distortions [17] and electron antineutrino back- 
ground flux |12j . It turns out that the constraint from 
deuterium photodissociation is the strongest one, so for 
this mass region we use the constraints on V-r. of [13j . 
where limits derived in |16j were updated. One should 
note that an idea of constraining of PBH production 
rate from an influence of PBH evaporations on the light 
element abundances had been proposed in pioneering 
paper [T8 . CMB distortion constraints obtained in [T7] 
are basically about the same order of magnitude, and 
limits from considering the neutrino background are 
somewhat weaker. For even smaller PBH masses, other 
useful limits exist, e.g., the constraints from hadron in- 
jection by PBHs during nucleosynthesis [H] (see the 
review [13), but we do not explore this region here. 

Numerical integration of ([3]) with the input given by 
© shows that for realistic peak with the finite width 
(particularly, E > 1), the good estimate for maximum 
value of Q.Qw(k) is 



QSwM « n GW (k ,r ) = 0.002 [ 9 -j L ) / ■ (V^) 2 



2E 2 



(9) 

(8) i.e., maximum value does not depend on E and is 
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approached at the scale ko (the particular shape of 
QQw(k), of course, depends on E - see analysis in the 
work [8]). Since the PBH constraints on weakly de- 
pend on E, formula @ can be used to obtain the limits 
on maximum possible Q,Gw{k) for second order GWs. 

Fig. [T] shows the resulting constraints on Qow(k) 
calculated using known PBH limits on the scalar spec- 
trum. The two curves correspond to the constraints 
obtained for two different gravitational collapse models 
(standard and critical) that assume different values of 
S c . Since for critical collapse larger d c was used, larger 
values of VLcwik) are possible in this case. 

It should be noted that in the case of a very narrow, 
(5-function like, spectrum of V-jz(k), the constraints on 

used by us, of course, do not hold anymore, and the 
formula ([9]) can not be applied. However, the analysis 
shows [7J that in this case maximal possible GW energy 
density is ~ 10 -7 — 10 -8 in the range of frequen- 

cies which we consider (~ mHz - kHz), i.e., the values 
are smaller than bounds of Fig. Q] So the presented 
bounds are the actual upper limits on flQw(k). 

Note that the constraints shown in Fig. [T]cover both 
the regions of / that are currently probed by ground- 
based experiments LIGO and Virgo (best sensitivity at 
/ ~ 100 Hz) and the ones that will be accessible to 
future space-based experiments LISA (/ ~ mHz) and 
BBO/DECIGO (/ - 0.1 Hz). 

3. THEORETICAL ESTIMATES FOR THE 
STOCHASTIC GWB 

Standard inflationary model. It has been derived 
in [19j that de Sitter expansion phase in the early uni- 
verse generates the stochastic GWB with scale invari- 
ant power spectrum through the mechanism |20j of am- 
plification of quantum vacuum fluctuations. So, in- 
flation is one of the most promising sources of GWB 
PHHH]. The prediction for it is given by PHHS] 

Here, tt r h 2 w 4 x 10 -5 is the total energy density 
fraction of radiation at present time, V;Jt is the en- 
ergy scale of inflation. From WMAP observations, 
V£} < 3.4 x 10 16 GeV, and fij^fc 2 < 2 x lO" 15 for 

/ > io- 16 Hz ng. 

GWB from preheating. The peak value of the back- 
ground density is estimated by the formula [571 HE] 

n P G W (k P eak)h 2 = firfc 2 -^- f^V /3 . (11) 
9 2 m P i \9*J 
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Figure 1. Constraints on the values of Qgw for induced 
GW background, calculated using formula (|9} and ex- 
isting PBH constraints on Vn(k). The forbidden val- 
ues of Qgw lie above the curves shown in the figure. 
Upper curve: constraints obtained using the assump- 
tion of critical collapse, with threshold value of density 
contrast S c — 0.45; lower curve: the case of standard 
collapse with S c = 1/3. 

It is assumed here that the inflationary potential at the 
end of inflation is given by 

i.e., /x is the effective mass [28] and the frequency of 
an oscillation of the field 4>, g is the coupling constant. 
Energy scale of inflation is V^f ~ ^fimpi, the Hub- 
ble function is H e — JVinf ~ M> an d the peak 
frequency (the predicted spectrum has a definite maxi- 
mum) is 



f peak « 6 x 10 1( y ^-flz ~ 6 x i ity J_Hz. (13) 

If fi is very small, fi ~ 10~ 18 mpi, then f pea k ~ 
200 Hz, i.e., the peak is situated in the region stud- 
ied by the LIGO detector. In such a case, however, 
the coupling constant g 2 is extremely small, of order of 
10 -30 . Note that in standard models of chaotic inflation 
the parameter \i is not free, u, ~ 10~ 6 mpi 1 and the peak 
frequency is ~ 10 8 Hz. 

The predictions of hybrid inflation models are more 
optimistic 1301 EH] • It was shown, in particular, that 
GWB can be relatively large in the LIGO/BBO region. 
For the case g 2 /X 1 the prediction of [30] is 

f peak « ^A 1 / 4 10 10 25 Hz, (14) 



4 



E. V. Bugaev, P. A. Klimai 



h 2 Q GW <10- 81 (^) . (15) 

Here, A and g are parameters of the inflationary poten- 
tial: A is the Higgs self-coupling, g is the Higgs-infiaton 
coupling. For instance, for A ~ 0.1, f pea k ^ 10 3 Hz, 
9 2 < 10 _ls one has h 2 il GW < 10~ 6 . 

In supersymmetric inflation models GWs can appear 
also as a result of an explosive decay of flat direction 
condensates [35]. Typical frequencies of these GWs are 
in Hz-kHz range, and the amplitude can be as high as 
h 2 n GW ~ 10" 6 [33J. 

First-order phase transitions in early Universe. 

The peak frequency is given by the formula [34; 

W „ 6 xlO-»n^(^) 1/6 (^)A. (16) 

Here, if* is the Hubble function at the time of GW pro- 
duction, /» is the characteristic frequency of GWs, T* 
is the temperature at the time of the phase transition. 
Typically, f,/H, ~ 10 2 [34], and f peak ~ 10~ 3 Hz if 
T* - 100 GeV and f peak ~ 10 2 Hz if T* - 10 7 GeV. 
The energy density at the peak frequency is [34] 

n% w { fpeak )h 2 « io- 5 (^) 1/3 (^) 2 , (17) 

i.e., ^q\y{ fpeak) ~ 10 9 . 

Cosmic strings. The estimate of the energy density 
is given by the formula [35] (in high frequency limit, 
/ > IO" 4 Hz) 




Here, G/i is the string tension, 7 is the radiation ef- 
ficiency, a is the initial loop size as a fraction of the 
Hubble radius. Pulsar timing data constrain ftg^(/) 
at / ~ 10~ 8 Hz, and, according to this constraint [36] , 
one has, at high frequency tail, 

Vow £ 10 " 8 - ( 19 ) 

Pre-big bang scenario. In string cosmology (which 
is the theoretical base of the pre-big bang (PBB) sce- 
nario [37]), the spectrum is growing with frequency, and 
the maximum value is determined by the value of the 
parameter M s /Mp (where M s is the fundamental string 
mass). According to the analysis of [38 , the maximum 
estimate of the energy density is 



This maximum value is reached at / ~ 100 GHz, i.e., 
in the region which is inaccessible for an experiment (at 
present time, at least). However, nonminimal variants 
of the scenario are theoretically possible [38 , in which 
the peak is shifted to lower frequencies without changing 
the peak amplitude (e.g., to / ~ 10 2 Hz). 

Brane-world models. In these models, in which our 
world is a brane embedded in a higher dimensional 
space, new geometrical degrees of freedom can be in- 
troduced (e.g., the position of the brane). Excitations 
of these degrees of freedom can lead, in principle, to 
detectable GW radiation. Naturally, the characteristic 
frequencies depend on the size of the extra dimensions 



GWB from PBH evaporations. GWs are neces- 
sarily produced during PBH evaporations in early uni- 
verse. It is shown in [40] [41] that the corresponding 
stochastic GWB can be substantial, depending on the 
parameters of PBH formation models. Unfortunately, 
typical frequencies of GWs are rather high (10 12 Hz or 
more) . 



4. CONCLUSIONS 

We have presented in the previous section that 
the typical predictions of the energy density of GWB 
(characterized by the value of £Igw) do not exceed 
10~ 8 — 10~ 9 , although the examples of nonminimal vari- 
ant of the PBB scenario and some variants of preheating 
models (see Eq. (|15[) ) show that larger intensities (10 -7 
or even 10 -6 ) are not, in general, excluded. 

The result, obtained in this paper, can be formu- 
lated in the following way: the contribution to the total 
GWB, which is induced by primordial density fluctua- 
tions, is constrained, in the region of ~ 10~ 3 — 10 3 Hz, 
as it is shown in Fig. [T] 

This bound should be taken into account in elabo- 
rating inflationary models because, in principle, the in- 
duced part of GWB (produced as a result of inflation) 
can significantly exceed the first order part. The well 
known example is the running mass model [H] . 

From the experimental point of view, the bound of 
Fig. [T] gives the maximum possible signal from second 
order GWs. It is important because if the experimen- 
talists will be lucky do detect GWs above this bound, 
it will mean that the signal is not due to induced GWB 
(e.g., nonminimal variants of the PBB scenario men- 
tioned above predict peak value of Slew close to 10 -6 ). 
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The bound shown in Fig. [T] is slightly more restric- 
tive than the nucleosynthesis bound [JH 03] (in LIGO 
region, at least). The latter bound depends on the effec- 
tive number of neutrino species N v while the former one 
is a straightforward consequence ot the PBH limit on 
the spectrum of primordial density perturbations (and 
Einstein's equations). 
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